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Abstract 

During the past decade, the interaction of high-intensity lasers with solid tai^gets has attracted much 
interest, regarding its potential in accelerating charged particles. In spite of tremendous progress in laser- 
plasma based acceleration, it is still not clear which particle beam quality will be accessible within the 
upcoming multi petawatt (1 PW = 10^^ W) laser generation. Here, we show with simulations based on 
the coupled relativistic equations of motion that protons stemming from laser-plasma processes can be 
efficiently post-accelerated using crossed laser beams focused to spot radii of a few laser wavelengths. 
We demonstrate that the crossed beams produce monoenergetic accelerated protons with kinetic energies 
> 200 MeV, small energy spreads 1%) and high densities as required for hadron cancer therapy. To our 
knowledge, this is the first scheme allowing for this important application based on an all-optical set-up. 
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Accelerating charged particles is of paramount importance in a wide variety of fields, rang- 
ing from medicine 1 1 and material science [|5D to being used to resolve the smallest structures 
of our universe dg, [ZD. The rapid development of high-intensity laser systems which are likely 
to exceed in the near future 10^^ W cm~^ [8, ^ rendered particle beam creation by laser-matter 
interaction feasible. Laser-driven accelerators offer the unique feature of ultra-high electric field 
gradients of several TV m~^, outperforming those in conventional accelerators by more than six 
orders of magnitude and thus offering the possibility of compact and low-cost devices ifioll . Elec- 
tron IjUlI] and proton beams have recently been generated by focusing high-intensity laser light 



2311 is real- 



onto solid targets. This mechanism of target normal sheath acceleration (TNSA) [[12- 
ized by the strong quasi- static electric field induced by the ionization and acceleration of electrons 
by the intense laser field. A further laser-plasma-interaction process, the skin-layer ponderomo- 
tive acceleration (S-LPA) resulting from the huge electric potential gradient of the plasma leads 
to ion beams of high density 124]. To date the energy resolution and the ion kinetic energy of 
the generated beams merely reach the parameters required for skin-deep neoplasm, i.e. excluding 
deep-seated tumors. Furthermore there is still controversy [|l5|l which beam quality, i.e. which 
total particle number, kinetic energy and energy spread will be accessible with the forthcoming 
multi petawatt-class laser systems M,^. 

In this article we investigate the vacuum post-acceleration of plasma-generated ions by means 
of lasers in a single or crossed-beams configuration. This scheme of laser particle acceleration 
was first proposed for electrons i25l- l27|l . The basic idea is to send the protons originating from a 
TNSA or S-LPA experiment through the crossing point of two laser beams at a half-angle 9 with 
respect to the z-axis (see Fig. lb for a scheme and a coordinate system). Within the simple plane- 
wave picture and assuming the laser fields to have the same amplitude, frequency, phase and to 
be polarized as shown in Fig. lb, the resultant electric field component along the symmetry axis 
of the set-up vanishes for all points on the x-axis. At the same time, the x-component, increased 
by constructive interference, violently accelerates the particles. Such a coherent combination of 
intense beams is experimentally feasible (see p. 52 of Ref. [8]). Subsequent motion of the charged 



particle ejected from the focal region and no longer interacting with the laser pulse may be taken as 
linear. The results of our theoretical simulations demonstrate that the protons gain kinetic energies 
larger than 200 MeV (employing two crossed beams each with a peak intensity of 1.9x 10^^ W 
cm~^) with an energy spread of roughly 1%. Furthermore (assuming a realistic repetition rate of 
10 Hz [18i|28|]), the total number of generated protons reaches 10^° min~^. For the first time, all 
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FIG. 1: Laser-proton acceleration schemes, (a) The protons, produced by laser-plasma interaction, are 
injected with the angle 9i with respect to the propagation direction of a pulsed beam through its focus. The 
laser field polarization is denoted by E and the propagation direction is given by k. The protons are ejected 
out of the focus in the polarization direction E. (b) Here, the protons are injected through the intersection 
point of two pulsed beams with crossing half-angle 9. The laser field polarizations are denoted by Ei and 
E2 and their propagation directions are given by ki and k2. The protons are ejected in direction of the 
resulting electric field E. 



requirements are fulfilled for broader radio-oncological use yj, |20 based on an optical accelerator. 



Results 



In order to generate ultra-strong accelerating fields |18|,|9D of 10^^ W cm^^one needs to focus the 
laser field to a beam waist radius on the order of the laser wavelength [|29n . This necessitates an 
accurate description of the fields beyond the widely-used paraxial approximation. The parameters 
of a linearly polarized Gaussian beam which propagates in the 2-direction and is polarized in the 
x-direction and subsequently rotated by 6i (Fig.[T^) will be used to model the fields, i.e. the beam 
waist radius wq, the Rayleigh length Zr = ttwq/X, where A is the laser's wavelength, and the 
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diffraction angle i?, e = wq/ Zr = X/ (vrwo)- The expressions giving the Cartesian field components 
E,.. E,. E^^ B^. B,. B.. a. well a. ,he expression for U,e power of Ire fields ,o order in^e 
diffraction angle can be found in the Methods section and for more details we refer to Refs. BOL 
3l|] . For the intensity profile of the employed Gaussian beams see Fig.|2l 

High-intensity laser systems provide their energy in short pulses which are already sufficient to 
accelerate particles to high velocities li32i1 . Employing pulsed fields also ensures that the particles 
injected into the focus get captured rather than reflected. To the lowest order in time, this can be 
described by multiplying the electromagnetic field components with a Gaussian temporal envelope 
factor, 

/ {t-z/cf\ ^ 
B ^ exp r-i^l (1) 



with At defined via the Full Width at Half Maximum (FWHM) pulse duration AtpwHM 



2-y/2 log 2 At. This approximation is valid for T/Atp^HM ^ 1? with T being the laser period. 
For the titanium-sapphire laser with wavelength A = 0.8 /im (T = 2.65 fs) and pulse durations of 
At > 10 fs used in our simulations, this turns out to be an adequate description. Hence, further 



temporal corrections [|33|] which describe the field solutions as a dual power series in the diffraction 
angle e and in the small ratio T / (211 At) can be neglected. 

The motion of an ensemble of N identical particles of mass m and charge q in the electric and 
magnetic fields E and B, respectively, of a laser beam is considered classically, with randomized 
initial distributions. The use of laser systems of high intensity (exceeding 10^^ W cm~^ for pro- 
tons) requires a relativistic treatment of particle motion. Thus, the dynamics is governed by the 
coupled Newton-Lorentz (or energy-momentum transfer) equations (given in SI units): 

^ = q {E{r,) + Ef- + x {B{r,) + Bf )) , 
^ = gc/3, ■ (Eir,) + Ef ) . 

(2) 

The relativistic energy and momentum of a given particle labeled with j are denoted here by 
£j = 'jjmc^ and pj = ^jmc(3j, respectively, with /3j = vj/c its velocity scaled by c, and 7^- = 
(1 — its Lorentz factor. The fields mediating inter- ionic interaction are given in the Methods 
section. To obtain the kinetic energy gained by interaction with a laser beam, numerical solutions 
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FIG. 2: Laser intensity and field strength, (a) Intensity profile of two crossed Gaussian beams in the 
propagation plane y = 0. For visualization purpose the crossing-half angle is set to = 35°. The brighter 
areas correspond to higher field intensity, (b) Vector field plot of the polarization plane z = at t = 0.2 fs 
for the single beam scheme (blue an^ows) with power P and crossed beam scheme (red aiTows) with power 
P/2 for each beam. The constructive interference of the crossed beams results in a higher electric field 
strength in the intersection volume. The background of the graph shows a density map of the electric field 
strength |E| of two crossed beams in the polarization plane. 

of the equations of motions will be sought. A numerical integration of equation ^ yields /3^" and, 
hence, 7^'^ at a later final time taken equal to many laser field cycles. Finally, one calculates 
the final kinetic energy of the particle from K^'^ = 

Foremost, we carry out simulations based on the coupled equations of motion equation ^ 
for an ensemble of 50 particles at later used particle densities in order to determine the dominant 
nearest-neighbor contribution of proton-proton interaction effects on the resulting particle beam. 
Due to the dominating ponderomotive laser forces which lead to a fast drifting apart of the en- 
semble's particles, for relativistic laser intensities, i.e. > 10^'' W/cm^ for protons, it turns out that 
the dynamics, the energy gain and its spread are influenced negligibly by inter-ionic interaction. 
This is in contrast to the non-relativistic laser regime. Here, the repulsive Coulomb interaction is 
the prevailing part of the interaction, which is in this case non-negligible compared to the electro- 
magnetic fields of the laser. As a consequence, the accelerated ions occupy a larger phase space 
volume. From radio-oncological point of view we are interested in proton beams of fully relativis- 
tic energies, hence for further calculations at these energies and densities it is sufficient to study 
the uncoupled equations of motions only, i.e. we set Bj'- = Ej'- = in equation (|2l). 

The definition of a coordinate system for the crossed beams set-up depicted in Fig.flh is given 
by the coordinate transformations xi = x cos 6* — 2 sin 9, yi = y, zi = xsinO + z cos 9 for the first 
beam and X2 = x cos 9 + zsin 9, y2 = y, Z2 = —x sin 6* + 2; cos 9 for the second beam, respectively. 
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The resulting field components which enter equation Q are 

= (El, + ^2,) cos e + {Eu - E2,) sin e, 

Ey = Ely + E2y, 

E, = i-Ei^ + E^a:) sine + {El, + E2,) cose, 
= {Bi,-B2,)sme, 

By = Bly + B2y, 

B, = {Bi, + B2.)cose. (3) 

Choosing a small crossing half-angle 6* leads to constructive addition of the dominating x- 
components of the electric fields in equation ([3]). The adding of further laser beams would further 
increase the laser intensity and hence the exit kinetic energy of the accelerated particles. With 
respect to the energy spread one cannot achieve substantial improvement, however. For all subse- 
quent simulations we restrict our analysis to the case of two crossed beams and set = 3°. 

In order to simulate a realistic particle injection into the focal point of the laser pulses, we con- 
sider an ensemble of 5000 particles initially randomly distributed in a micron-scale volume Vf^^, 
oriented along the 2-axis. Initially, we have to ensure that the ensemble is not already exposed 
to the laser fields. This is realized by starting the simulations at time t < —5 At. At these initial 
times the electromagnetic fields are damped by the Gaussian envelope factor such that the parti- 
cles' motion is only negligibly influenced by the external laser fields. Fig. 3 compares the motion 
of a proton starting at the same spatial point but having two different initial times. One can see 

1.5 /' ~\ 
^1.0 / 

^ / \ 

0.5 / a\ /''\ ^ 




FIG. 3: Proton dynamics. Kinetic energy K for one proton initially being at rest and located at x = z = 
A/30 and y = 0. The initial time is = (dashed line) and ti = —10 At (full line). The single proton 
dynamics is compared at same laser parameters. 
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that in the unphysical case when the particle is directly exposed to the laser fields (dashed line), 
it immediately gains energy and is ejected out of the focus after one laser cycle. In a realistic 
setting, we choose an initial time oft = — lOAt: the particle slowly starts to oscillate and then 
gets captured and accelerated by the approaching pulse (full line). Choosing the initial time to be 
t = would lead to an energy gain overestimation by a factor of approximately three. 

The volume initially containing the particle ensemble has a length of the order of the laser 
wavelength and a radius of tens of nanometers, dependent on the focus diameter of the applied 
laser system, ensuring that all protons are exposed to a homogeneous field. The particles will be 

TABLE I: Variation of laser parameters. Average particle kinetic energy K and its percentual spread for 
different laser system parameters. Ni = rii • Ffo^.^, is the number of ions one can accelerate as one bunch 
with n^^''* 10^^ cm~'^ and n™^^ k, 10^^ cm~^ is the ion density of the source used and Vf^^^ denotes 
the volume initially containing all ions. The crossing half-angle is ^ = 3°. The optimal particle injection 
angle for the single beam set-up is 9i = 3° for the S-LPA source and 9i = 50° in case of the TNSA source, 
respectively. For two crossed beams the paiticles are injected with an angle 9c with respect to the symmetry 
axis (z-axis) of the laser beam configuration. In case of the S-LPA source we have 9c = 0° and for the 
TNSA source 9c = 50°. 
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28.1 ± 1.2% 
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crossed 


59.4 ± 1.0% 
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S-LPA source, P=40 PW, At= 


■lO.l fs, Wq = 
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single 


113.2 ± 1.6 % 
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crossed 


233 ± 1.0 % 






1.0- lO'' 




S-LPA source, P=100 PW, At: 


=23.8 fs, wq = 


: 2A 




single 


73.2 ± 1.6 % 






1.3-10^ 


crossed 


152 ± 1.0 % 
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TNSA source, P=100 PW, At= 


= 14.4 fs, Wq = 


: 2A 




single 


64.6 ± 0.7 % 






I.O-IO^ 


crossed 


141 ± 0.5 % 






1.0-10^ 
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assumed to possess initial kinetic energies distributed normally around a mean value K and having 
a spread A/T . As a source we take protons originating from laser-plasma interactions, such as the 



S-LPA mechanism, with K = 11 keV ll24ll and assuming a large energ y sp read of = 100 % 



or from the TNSA mechanism, with K = 1.2 MeV and AK = 25 % [ 



13h . The total power and 



the beam pulse duration are varied. Note that the peak intensity of one linearly polarized 10 PW 
laser beam focused to wq = IX is already / ~ 9.6 x 10^^ W cm~^ [8, ^. 



Discussion 

Our main interest is in the energy gain, or exit kinetic energy, of the nuclei, their trajectories 
and, hence, the aspects that determine the quality of an accelerated beam of such nuclei. In Tab. 1, 
simulation results for the laser acceleration are summarized. The single and crossed beams scheme 
are compared at same total laser power P, pulse duration At and focus radius Wq. The energy gain 
is ranging from 59 MeV to 233 MeV in case of the crossed beams setup and from 28 MeV to 1 13 
MeV for injection of the ensemble into the focus of only one beam, respectively. Looking at the 
energy spread one can see that for the crossed beams it is always 1% and a little higher for the case 
of the single beam setup. 

Moreover, one can see from Tab. 1 that for the systems studied here, the average exit kinetic 
energy K obeys for constant injection angle and injection energy the rough scaling behavior K oc 
/ oc P/wq. This behavior results from the fact that the optimal acceleration regime depends 
strongly on the pulse duration At rather than only on the electric field strength. It is achieved for 



the laser-particle interaction length being of the same order of the Rayleigh length. In Ref. [|34l] 
the same scaling behavior was derived for electrons. Therefore, in order to maximize the energy 
gain of the accelerated protons we first choose the laser power P and the focus radius wq, and then 
adjust to the optimal pulse duration At. 

A further scaling law can be derived for the particle number, which is proportional to the focal 
volume. Using Vt^^,,, (x. ■ Zr and the definition of the Rayleigh length Zr oc Wq, one obtains 
Kocus Wq. The typical particle number needed for ion cancer treatment is 10^-10^° per shot with 



a repetition rate on the order of 5 Hz 113511 . depending on the ionic species. For a typical interaction 
volume, one needs at least an ion density of rij = 10^/l^o,„, ~ 10^° cm^'' in order to render our 
scheme feasible for medical applications. Using plasma-generated protons one obtains with the 



TNSA mechanism an ion density up to the order of n™^^ ~ 10^^ cm ^ and the S-LPA ion source 
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generates a density of 



10^^ cm ^ n24\ . The latter yields particle numbers of 10^ per laser 



2811 this ion number is sufficient 



shot. Combined with lasers operated at 10 Hz repetition rate [81, 
for cancer therapy while for the TNSA mechanism a modest improvement would still be necessary. 

Such improvement can e.g. be achieved by substituting the assumed titanium- sapphire laser 
by a super-intense CO2 laser with a typical wave length of A = 10.6 /im. Using the fact that 
the waist radius wq oc A and hence the Rayleigh length Zj. oc A, the focal region thus increases 
by three orders of magnitude. Consequently, the ion density needed decreases by three orders of 
magnitude and the needed laser intensity by two orders of magnitude. In the near future the use of 
high repetition rate laser systems M, \2M will further decrease the number of protons needed per 
shot. 

The combination of the laser-plasma mechanism as a proton source and the post acceleration 
process by means of single or crossed laser beams places laser acceleration of particles on the 



cusp of medical feasibility utilizing present-day or near-future laser technology 



m 



anticipated costs are presently on the scale of those for current synchrotron facilities (cf. [|36|l and 



281, while 



13711 ). The rapid advancement of laser technology renders significant reduction likely for the near 
future. All requirements needed for broader radio-oncological use may be achieved: sufficient 
proton density and exit kinetic energies, and sharp energy spread of approximately 1 %. The 
scheme that we introduce in the present work calls for tight focusing mechanisms n2% and relies 



on results of laser-plasma- interaction research 1112 
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Methods 



Description of a tightly focused linearly polarized beam 



The electromagnetical field components of the tightly focused linearly polarized beams which we 
employ in our simulations can be given as a power series in the diffraction angle e. Their functional 
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form is given by 
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With uj the angular frequency of the fields, ^ = x/wq, v = y/wo, C = z/zr, iI}g = tan ^ C, and 



r = \J + p = r/wo, the remaining symbols in equation dH) have the following definitions 



E 



where 



(—) cosiijj + nipc); n = 1,2,3,- ■■ 
\ w / 

Sn = (—) sm{ip + niljG), 
\ w / 



ip = ipo + ut - kz - —; R = z + —, 
2R z 



(5) 



(6) 



and ^0 is a constant initial phase. Also, t is the time and k = 27r/A is the wavenumber. On the 
other hand, with Eq Eqi, the power expression may be given, to the same order in e as the field 
components, by 



Pi 



4 c/io 
45 (e 

+y(2 



i) +2(5) +6(2, 



195 /e\io- 
~ V2 



(7) 



where c is the speed of light in vacuum, /iq is the permeability of free space and Eqi is the electric 
field amplitude, with / standing for linearly polarized. Note that Eqi oc \/Pi and that the leading 



term in Eqi is inversely proportional to Wq. For more details we refer to Refs. 
Relativistic coupled equations of motion 



31II 



The fields mediating the inter-ionic interaction in the relativistic coupled equations of motion 
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@ are modeled by Ef = Y.k^j {-^(f^jk - | ^jfe) and Bf = Y.k^j (V x A^k) with j, k e 
{1,2,...,A^}. The interaction potentials read 

9 1 



47reo \rj - n 



(8) 



A,. = -^f.. + ^>^f^V (9) 

with the relative displacement Vjk = rj — and eo being the vacuum permittivity. Equation dS]) 
is the scalar part of the interaction given by the Coulomb potential, whereas relativistic effects 
such as retardation and current-current interaction are included in the Darwin vector potential up 



to 113 8h in equation Typical kinetic energies of the accelerated protons are of about 200 

MeV (cf. Tab. 1), which corresponds to ~ 0.3. Consequently, the truncation of the interaction 
up to is justified and higher-order contributions such as those up to in Ref. [,39J will 

not be taken into account. 

Equation dll) is a differential algebraic equation of the form 

/i ri, . . . , T-AT, /?!,..., Pn, —rr^ —rr: i 



dt V dt' ' dt 



= /jv I • • • , T-JV, Pi, • • • , Pjv, • • • , t 1 . (10) 

For each time step the system has to be solved algebraicly with respect to . . . , Then, a 
standard fourth-order Runge-Kutta algorithm is used to integrate the remaining system of ordinary 
differential equations. 
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